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ABSTRACT 

Under the separability assumption on the augmented density, a distribution function can be always constructed 
for a spherical population with the specified density and anisotropy profile. Then, a question arises, under what 
conditions the distribution constructed as such is non-negative everywhere in the entire accessible subvolume 
of the phase-space. We rediscover necessary conditions on the augmented density expressed with fractional 
calculus. The condition on the radius part R(r 2 ) - whose logarithmic derivative is the anisotropy parameter - is 
equivalent to w~ l R(w~ l ) being a completely monotonic function whereas the condition on the potential part is 
stated as its derivative up to the order not greater than | -/?n being non-negative (where /3q is the central limiting 
value for the anisotropy parameter). We also derive the set of sufficient conditions on the separable augmented 
density for the non-ne gativity of the distribution , which generalizes the condition derived for the generalized 
ICuddefordl system by ICiotti & Morgantj (l2010h to arbitrary separable system s. This is applied for the case 
when the anisotropy is parameterized by a monotonic function of the radius of Baes & Van Hesel d2()07). The 
resulting criteria are found based on the complete monotonicity of generalized Mittag-Leffler functions. 



1 . MODELS FOR SPHERICAL DYNAMICAL SYSTEM 
1.1. Distribution function 
Suppose that 'Fir, v\t) is a phase-space distribution so that 



X 



is the number of tracers in any measurable phase-space vol- 
ume S at time t. Here r is the position vector in the configu- 
ration space and v = r is the velocity. We only consider the 
system in equilibrium and thus the distribution must be time- 
independent. The distribution of a spherically symmetric pop- 
ulation in a steady state is also invariant under transforms in 
SO(3) so that T{r;v\t) = T(r;v r ,\\v t \\) where r = ||r|| and 
t = r/r are the radial distance and unit vector while v r = v • r 
and v t = v - v r r are the radial and tangential velocities. If we 
adopt the canonical spherical polar coordinate (r, 6, (p), they 
are given by 



(v 2 + v 2 )>; 



where (v r , vg, tu) = (r, ri), r<f> sin 6) are the velocity compo- 
nents projected onto the associated orthonormal basis. 

In order for the distribution to be indeed time-independent, 
it must be invariant under dynamic evolutions of tracers, that 
is, the distribution is a time-independent solution to the Boltz- 
manrQ transport equation. For typical stellar dynamical appli- 
cations, the trajectory of each tracer is its orbit under the exter- 
nal potential, which may or may not be self-consistently gen- 
erated by the tracer population. The transport equation for this 
case results in the collisionless Boltzmann equation (CBE), 
whose solution is c omple tely characterized by the theorem 
due to J. Jeanfl The I Jeansl theorem indicates that if the given 
time-independent spherically-symmetric distribution function 
(df) is a solution to CBE with a generic static spherical poten- 
tial O(r), it must be in the form of 

T = T{S,L 2 ) 

1 Ludwig Eduard Boltzmann (1844-1906) 

2 Sir James Hopwood Jeans (1877-1946) 



where 



£ = T'(r) - \v 2 



L = ||L|| = rv t , 



are the two isotropic isolating integrals admitted by all such 
potentials, namely, the specific binding energy and the mag- 
nitude of the specific angular momentum, respectively. Here, 



$Oout) - <t>0) 

(J>(oo) - O(r) 
-<P(r) 



if r out is finite 

if '"out = 00 and |0(oo)| < oo 
if /"out — 00 and O(oo) — > oo 



is the relative potential with respect to the boundary r out . The 
system not bounded by a finite boundary radius is represented 
by '"out = 00 with <J>(oo) = lim^K, <£>(r). If r out or (D(oo) is 
finite, then T{& < 0, L 2 ) = because by definition S > for 
all tracers bound to the system (and bounded by r < r out ). 

1.2. Augmented densities of a spherical system 

Integrating the spherical two-integral df T(&,L 2 ) over the 
velocity space results in a bivariate function of and r 2 , 



NOP, r 2 ) = jjj d% T(S = ¥ - \v 2 , L 2 



r 2 v 2 ), (1) 



which is referred to as the augmented density (AD). The inte- 
gral here is formally over the whole velocity subspace, but 
if '"out or <D(oo) is finite, it is essentially within the sphere 
v 2 < 2 X P since <F(S < 0, L 2 ) = for these cases. With »P(r) 
specified, the AD yields the local density v(r) via 

y(r) = N[*P(r), r 2 ]. 
Similarly, the augmented moment functions are given by 



m k , n (^,r 2 ) 



T(& = x ¥- jV 2 ,L 2 = r 2 v 2 ) 



in Jj iv, do, W'lFfT - A?). (2a) 



([i 2 <2Y) 

!V>0,Dt>0 
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Changing the integration variables to (£, L 2 ), these are repre- 
sented to be a set of integral transformations of the df, 



2n 

r 2w+2 



2n+2 IIr 

II 



dSdL 2r K k -?L 2 "T(S,L z ) 



dS dL 2 0(7C) \K\ k -s L 2n T(S, L 2 ). (2b) 



'' " J.h: ' 2 >0 
Here 0(x) is the HeavisideQ unit-step function and 
fO if r out or <l)(oo) is finite 



£0 = 



-00 if lim^oo *P(r) = -<J)(oo) — > -00 



is the lower bound of the binding energy. The transform ker- 
nel and the domain in (&, L 2 ) space over which the integral is 
performed are given by 

<K(&, L 2 ; x ¥, r 2 ) = 2QV - &) - r 2 L 2 , 
T = { (S,L 2 ) I £ > & lh L 2 > 0,7C > 0} . 

Note 7C is u 2 expressed as a function of 4-tuple (£, L 2 ; *P, r 2 ). 

2. MATHEMATICAL PRELIMINARY 

2.1. Fractional calculus 

Definition 2.1 77ze Riemann^-Liouvill^ integral operator of 
arbitrary non-negative real order A > is given by 

[m » = ») 

where T(x) is the gamma function. 

This is a trivial generalization of the Cauch>0 formula for re- 
peated integrations. For < A < 1, this is also recognized as 
the generalized AbeQ transform with the classical case corre- 
sponding to the A = I case. We also define 



Definition 2.2 the fractional derivative for A > such that 

f(y)dy 



x J dx^ A ^ aJx 
1 



7 
d ui 



r(l-{A}) dx m 



J a (X- 



dx A 



f (A) (x) 



(0 < \A) < 1) 



(Ml = 0) 



(4) 



where \A~], [A], and {A} — A — [A] are the integer ceiling, the 
integer floor and the fractional part of A, respectively. 

Note equation is a generalization of the differentiation for 
positive real order as is equation @ of the integration. These 
definitions extend to include a negative index using 



Definition 2.3 for arbitrary real A 

\I A f — + a d x A f and vice versa. 



3 Oliver Heaviside (1850-1925) 

4 Georg Friedrich Bernhard Riemann (1826-1866) 

5 loseph Lioville (1809-1882) 

6 Augustin-Louis Cauchy (1789-1857) 

7 Niels Henrik Abel (1802-1829) 



The basic composite rule for the Riemann-Liouville opera- 
tors is that, for any pair of non-negative reals A and £ 



. by i 

210 i 



(6) 



which may be shown by direct calculations using the Fubinfl 
theorem and the EuleiQ integral of the first kind for the beta 
function, that is, 



["dyix-yf 1 P 

J a J a 

= r dwf(w) f 
-I' 



dyix-yf' 1 I dw{y-w) A - x f{w) 



dy{x-yf l (y-w) 



dw f(w) (x — w)' 



Next for any real A and a non-negative integer n 



fdrCi-^-V" 1 . 
Jo 



±+f* f= +f i - i f. il+f A f = + f A -"f (7) 

feaJx J aj x J > A v n a Jx J a J x J y ' 



dx" 



Here the latter follows the former (« = 1) by means of in- 
duction. The n — 1 case is proven by direct differentiation of 
equation ([3]i for A > 1 and the fundamental theorem of cal- 
culus for A = 1 while the same case with A < 1 is essentially 
trivial from the definitions of fractional derivatives in equa- 
tions (01 and ©. Together they also indicate that 



(8) 



for non-negative reals A, g > and arbitrary function fix), 
provided that all the integrals in their respective definitions 
absolutely converge. Next we observe for A > that 



+ r A+1 f = +C A f ( x ~ fl )' 1 /( fl ) 
a Jx ■' ~ "Jx ■> ' 



(9a) 



iX^ + i) 

thanks to the fundamental theorem of calculus (A = 0) and 
integration by part. By means of induction, this generalizes to 



(x - a) A+k f (k \a) 



+ C A+n An) _ + r A r _ y (x-a) j 

aj x J aj x J 2-1 Y(A + k + 1) 



(9b) 



where n is any non-negative integer, and we also find that 

d " +C A * _ + C A An) , V (* - a)f- k f- k \a) 



dx' 



+ C* f = + C*f(n) y Kx-ap 

n aJx J "Jx J Z_i Y{\ + 



r(i + A - k) 



(10) 



for A > and any non-negative integer n. The last implies that 
fractional derivatives in equation © are alternatively given by 

+y f _ + , m -A in) (x-at- A f k \a) 



d*W- 



r(i + k - a) 



where A > and n = 0, 1, . . . , \A]. 

Using these and equation ( TTOb . we can also derive that 



UJ 



(5) tf x U A A f = ://"7 - 2 q, :d/- k f(a) (x - a} 



k=\ 

UJ 



^:d x A f=tdt A f-Y J c i 



k=\ 



_ + a d x A - K f(a) 
(x-a) k+ f 



(12) 



8 Guido Fubini (1879-1943) 

9 Leonhard Euler (1707-1783) 
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for non-negative reals A, £ > and arbitrary function /(x), 
provided again that all the integrals in their respective defini- 
tions absolutely converge. Here Ci k are given by 

am 



c* - l - 



r(i +0 

(-l)^ J+ *((5) 



LfJ+* 



(+ case) 



(- case) 



m -5) 

where < 5 = £ - |_£J < 1 is the fractional part of and 

(aX = U%i(a + 1 - j) ; Kt = n "=i(« - 1 + J) 

are the Pochhamme£3 symbol. The falling product (a)~ fol- 
lows the combinatorist's convention whereas the rising one 
(a)+ does the analyst's. Note these are related to each other, 

(-a)- = (-1)»„ + ; (a)~ = {a-n+ I)* 

and also to the gamma functions, 

F(a + n) _ F(l + a) 



(a)Z = 



rxa) 



r(l+a-n)' 



The last may be used to generalize the Pochhammer symbol 
for non-integer n. Together equations © and ( fT2b provide the 
generalization of equation © for any pair of reals £ and A. 
The simplest specific result of fractional calculus would be 

Lemma 2.4 for real A and a > 0, 

T(a) x a+A - x 



F(a + A) 

This is formally a generalization of the result, namely 

d"x a 



dx" 



= (ay„x a 



(71=0,1,...) 



(13) 



(14) 



although the last is in fact valid for any a. 
We formalize an obvious but important fact, namely 

Lemma 2.5 for A > and x > a, if f > in [a,x], then 

af x f( x ) - 0- Moreover *f f + provided that the support 
of f in (a, x) has non-zero measure. 

Next, if a is finite, then for £ > 

+ rf, (*-a) f 

while for < A 

+ af!f = + ad'f 



f dff^7[x-(x- fl )f], 

1 (s) Jo 

£ < 1 and n = 1 , equation dTTb results in 

/(«) 



rxi-A) 



(x - a)"- 1 + 



It then follows that 
Lemma 2.6 /or a + +oo, 



lim 



#//« f(a) 



(x-ctf r(^+i)' 



(15) 



which is valid for £ > if f(x) is right-continuous at x — a 
or for £ > -1 if /(x) is right-differentiable at x = a. Equa- 
tion (15[ for £ = is equivalent to the definition of the right- 
continuity while for £ = -1, it becomes lim t _, +(x — a)f (x) = 
which holds if /'(a) is finite. Equation $15[ implies that 

10 Leo August Pochhammer (1841-1920) 



Corollary 2.7 iff(x) is right-continuous at x — a and f(a) is 
finite, then f(a) = 0/or A > 0. 

Next we examine the behaviors of fractional calculus oper- 
ators under the Laplac^3 transform. For this, we first note a 
general property of the Laplace transform of the derivative, 



s" +1 £ [f(x)] = L [f" +l> (x)] + £» sJf"-H0), (16) 

X— IS X— is 

which is valid given that lim x _,oo G~ lx f^ n \x) = for suffi- 
ciently large s (which is required for the Laplace transform 
to converge). Equation ( fT6l ) is proven via integration by part, 

f dxe~ sx ^^ = -/(0) + s [ dxe~ sx f(x) 
Jo dx Jo 

for n — and the induction completes its proof for any non- 
negative integer. In order to generalize equation ( TToT ) to in- 
clude the fractional derivative, we next consider for A > 

dxe" dy{x-y) A - l f{y) 
o Jo 

dyf(y)\ dx(x-y) A - l e-" 

Jij 

dyf(y)—r duu 

i Jo 

With the Euler integral of the second kind for the gamma func- 
tion, we find that 



X— 1 su 



s A £ \tf x A f(x)] = -C [fix)]. 



(17) 



The Laplace transform of an arbitrary real-order derivative is 
then found by combining equations ( fT6l ) and dTTb . 

2.2. Post-Widder formula 

Theorem 2.8 (Post-Widder) If 4>{t) is continuous for t > 
and there exist real 3 A > and 3 b > such that 

e~ bt \^)(t)\ < A for v f > 0, 

then the Laplace transform 



fix) = £ [0(f)] 



•f 

Jo 



dfe^V(f). 



(18) 



converges and is infinitely differentiable for x > b. Moreover, 
<p(t)for t > may be inverted from f(x) usi ng the differential 
inversion formula SP^sTJ^WddMrnlll . 

i (—1)" /7I\«+1 i„,/fl\ 

0(0 = -C' 1 [fix)] = lim L-i. (-) /">(-). (19) 

In literature, the last formula is typically named after E. PosFl 
or together with D. WiddeQ A rigorous proof, which is be- 
yond the scope of this paper, may be found in a standard text 
on the Laplace transform. However its heuristic justifications 
abound and are easy to observe. For instance, direct calcula- 
tions using equation ( TT8l indicate that 



f"\x) = i-irf^dtfe-^it) = J"d**»e-V$. 
ii 



Pierre-Simon Laplace (1749-1827) 

12 Emil Leon Post (1897-1954) 

13 David Vernon Widder (1898-1990) 
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and thus we find that 
(- 



Note then 



where 



P(s; n) = — e 
n! 



is the probability density of the PoissorQ distribution with 
a mean of s = n. It follows that as n — > oo, the relative 
dispersion decreases and so (p(st/n) — > <p(st/n) = 0(f), which 
results in the Post-Widder formula. Note however that the 
convergence of equation (TT9b by itself does not necessarily 
imply that f(x) is the Laplace transformation of <p(t), which is 
rather a part of the condition for the formula to be valid. 

2.3. Completely monotonic functions 

Definition 2.9 A smooth function f(x) of x > is said to be 
completely monotonic (cm henceforth) if and only if 

(-l)7 (n) (x) > (x>0,n = 0,1,2,...)- (20) 

The definition extends to x > if fix) is right-continuous at 
x = 0. Some basic properties of cm functions are: 

Lemma 2.10 Let f and g be cm. Then, 

1. (-1 )"/™ for any non-negative integer n is cm. 

2. If F > in (0, oo) and f — —F', then F is cm. 

3. J x f(y)dy is a cm function of x if it converges. 

4. af +bg where a and b are non-negative constants is cm. 

5. f ■ g is cm. 

6. If F > in (0, oo) and f — F' , then g o F is cm. 

7. exp(/) is cm. 



Items I and 2 are essentially trivial from Defintion 12.91 and 
item 3 is simply a particular case of item 2. Item 4 follows the 
linearity of differentiations while item 5 is shown using the 
LeibnizT^rule, that is, (here KM is the binomial coefficient) 

<ntA = yM H) ^ Hr ^ (21) 

1 L> dx" Mk ( ' dx k ( ' d*"-* - ( ' 

k=i) 

The last two may be shown using the Faa di BruncH formula 
(i.e., the generalized chain rule), 

(g o F) w (0 = Zlo9 (k) [F(t)] ■ B„, k [f(t)J'(t), ■ - ■ ,/ ( "^(0] 

(22) 

where F'(t) = f(t) and B n # is the BelFl polynomial, that is, 

/> *" i,v - ■ • ■ ' v " a) = ( r ( it) ;o ( 2t) ;i " ■ • 

Here the summation is over all sequences (jo, ji, . . .) of «on- 
negative integers constrained such that 

Em=0 7m = *: ; Zm=o('« + Dim = »• 

14 Simeon Denis Poisson (1781-1840) 

15 Gottfried Wilhelm Leibniz (1646-1716) 

16 Francesco Faa di Bruno (1825-1888) 

17 Eric Temple Bell (1883-1960) 



2m=0 "Urn = U - k 



and thus j m > indicates that j m = for m > n - k (n.b., if 
otherwise, 2m=o"Vm > n — k, which is contradictory). The 
property 6 follows this because 

n-k - 2 m= jlm+l = 2 2 m=0 m ( j 2m + J2m+l) 

is even. That is to say, if / is cm, the parity of the Bell polyno- 
mial in equation ( 1221 is (— l) n ~ k , and thus, given that g is also 
cm, the parity of every term in the sum on the right-hand side 
of equation ( 1221 ) is (-1)". Equation ( 1221 also indicates that 



d"exp[/(Q] 
dt" 



exp[/(f)] ■B n [f'(t)J"(i), 



,f"- k+1 \t)] 

(23) 

where B„ is the «-th complete Bell polynomial, that is, 
B n (x\, . . . , x n ) = Y!k=\ B n ,k{xo, . . . , x n -k). 

Note 

" ~ Dm=0 72m = 2 Em=0 m (72m-l + 72m)- 

is even. Hence if / is cm, the parity of the complete Bell 
polynomial in equation (|23| ) is (-1)" and so exp(/) is cm. 

The archetypal example of a cm function is fix) = e~ x . 
Other elementary examples of cm functions include: 

1. /(f) = r 6 (f > 0) is cm if and only if 6 > 0. 

2. /(f) = ln(l +T 1 ) is cm. 
These are proven through 



d"x 



= (-S)- X ~ s ~ n = (-1)"^- 

dx n ( )n 1 j x"+ s ' 



d" +l ln(l +x"') 
dx n+l 



(-l)" +1 n! 



1 

v-n+l 



1 



(1+x) 



n+1 



(24a) 
(24b) 



Following this and Lemma l2.10l are 

Corollary 2.11 Let g(t) be cm, then both t~ s g{t) with 5 > 
and g(t p ) with < p < 1 are cm. 

proof. The first is obvious thanks to Lemma \2.10\ 5. The last 
follows Lemma \2.101 6 with F(t) = t p since F' = pt p ~ l for 
< p < 1 is cm. q.e.d. 

Corollary 2.12 For < p < I and a, b > 0, these are cm: 

fit) = r a (i + t p y h ; /(f) = r fl (i + r p ) b . 

proof. Let F(t) = c + t p . Then F' = pt p ~ l is cm for < p < 1. 
Hence first (g o F)(f) = (1 + f'') - * with c = 1 and gr(iy) = ur fe 
/or < p < 1 b > is cm. Next, with c — ant/ 9(10) = 
Mn(l + w~ l ), we find that (g o F )(f) = Mn(l + r p ) is cm for 
< p < 1 anc/ > 0, a«d so is (1 + r p ) fc = e xp[fc ln(l + r p )]. 
The final conclusion follows Corollarv \2.1 71 q.e.d. 

The fundamental result characterizing cm functions 
dBernsteinll 1928k IWidder|[l94Tl) is due to S. BernsteirQ 

Theorem 2.13 (HausdorfF-BcrnsieiTi-Widder) 

A smooth function f(x) of x > is completely monotonic if 
and only if f(x) — L oT xt dp(t) where p(t) is the Borel mea- 
sure on [0, 00 ), that is, there exists a non-negative distribution 
4>(t) > oft > such that equation ( 1781 ) holds. 

18 Ceprefi HaTaHOBHi BepHinTCHH (1880-1968) 
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The 'if'-partis elementary since 



f n \x) = (-1)" dtfe-^t) = (-1)" £ [?'>(?)]. 

Jo »-»* 

Although the complete proof of the 'only if -part is beyond 
our scope, the partial proof follows the Post-Widder formula. 
That is, if the inverse Laplace transform (f>(t) = -C^,[/(x)] 
of a cm function fix) is well-defined, then equation ([T9T i, 
provided that it converges, indicates that <p(t) must be non- 
negative in the positive real domain. 

2.4. Miscellaneous 

We note an additional auxiliary relation, which will be used 
throughout this paper: that is, for any non-negative integer n 
and arbitrary differentiable function f(x), 



( x2 iJ (xf) = 



x" 



+1 d"(x"/) 
dx" ' 



(25) 



which may be proven via the induction on n fsee lAnlHoTTbl 
theorem A3). In fact this is also equivalent to a lemma 



*%+!)(*) = ^-[x" + 7(„)(x)j 

ax 



(26) 



where 



f(n)(x) = 



dx" 



This lemma may be proven directly via 



f(n+l)(x) 



d" 
dx" 
d"(x"/) 
dx" 

1 + n 



d(x ■ x n f) 



dx 



d" 


x"f + X 


d(x»/)l 


~ dx" 


dx 


ln\ d k 


x d"- k 


rd(x"/)i 


\kj dx* dx"-* 


dx 



( 1+ " + X dx) 



d \ d"[x"f(x)] 1 d 



dx" 



x" dx 



/(«)(*)] 



where we also used that d*x/dx* = if k > 2 and the Leib- 
niz rule (eg. [2Tb. The theorem in equation ( T25] l implying the 
lemma in equation d26T ) has been shown in I An! d201 lbl corol- 
lary A4) whereas the opposite implication may be deduced 
because the induction step for the proof of equation ( T25T) fol- 
lows equation d26l ) as 



d" +1 (x" +1 /) 
dx" +1 



1 d 

x" dx 



, , d \n I 1 / , d \" +1 



Fractional calculus also generalizes the lemma in equation 
generalizes. In particular, for a non-negative integer n 
andO < 6 < 1, 



tC\x n+S f) = 



r(i 



l fV 



f(y) dy 



yf 



jcrv + 7) 



x« +1 r 

~ m -5) J 

i r l dtf 
r(i-tf)J ~ 



1 t" +d f{xt) dt 
(1 - t) s 

n + S d" +1 [x" +1 /(xf)] 



(27a) 



dx' 



•n+1 



— f 

(1-5) Jo 



Jc i/" + % +1) ( J /)di/ 



x« +1 r(i 

x n+l OJx 



(x - y) s 



x" + % + i)(x)] 



x n+d /(„ + i)(^) = c^"V + W + V + 7)] 



(27b) 



(27c) 



Note that for the 5-0, the last results in equation ( 1261 1. The 
middle for the same case is consistent with the fundamental 
theorem of calculus given equation d26l i indicating 

-*" + 7(«)M = *^/kWt^ + f ffwM d y (28) 

provided that /(„)(0) is finite. Equations ( 1261 ) and d27b imply 

Corollary 2.14 /or a non-negative integer n, if /(„+i)(x) > 
/or x > 0, men tdx'ixff) > 0/or x > ana! n < fx < n + 1. 

In fact, the successive applications of this with a descending 
subscript furthermore suggest that, if /„)(x) > for x > 
and a non-negative integer n, it follows that idj'ixff) > for 
x > and a ny v p < «. 

Corollary 12.141 with an integer p may be generalized alter- 
natively, namely, 

Theorem 2.15 for a non-negative integer n, if x a /(„ + i)(x) is 
cm, then x fl /„)(x) is also cm. 

proof. If x a /„ + i) is cm, men by the Bernstein theorem, there 
exists a non-negative function h(u) > of u > such that 



r 

X a fn+\){x) = I 

Jo 



due~ x "h(u). 



The complete monotonicity of x a fi n ) can then be shown di- 
rectly using equation ( I2.SI ), which indicates that 

x fl / ( ,o = x"-"- 1 dyy n f( n+1) (y)= dtt" a \ duer m m, 
Jo Jo Jo 



dV/ W ] _ , 
dx* V ^ 

Finally, we also note 



f df?" +/; - fl r 

Jo Jo 



du e- x,u u k h(u) 



Lemma 2.16 for a non-negative integer n, iff^" +l \a) is finite 
and f (0) (a) = ■■■ = f k \a) = 0, then + a d" +5 f(a) = 0/or 
0<6<\. 

proof. Here we assume a — 0, but the similar argument holds 
for any finite "a " accompanied by a simple translation. First, 



°j* J m-s)j 



f(xt)dt 
(1 - tf ' 



4- <i n+6 r 

/ : 



1 



o 



l d n+ W s f(y)] 



dy" 



•dt 



Here the latter follows the former because 

d" +x [x l - s f{xt)] = tn+s d" +i [y'- s f{y)] 
dx" +1 

Finally, given the Leibniz rule, 
d B+ V _ 7(«/)] 



(29a) 
(29b) 



dy" 



y=xl 



= y x - s f (n+l \y) 



dy" 



/ w (y) 

-* ,,n+S-k ' 



which identically vanishes for y — if the condition part of 
Lemma \2.16\ with a — holds. Here the conclusion follows as 
the integrand of equation \29b\ with x — is also zero, q.e.d. 
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3. FRACTIONAL CALCULUS ON THE AUGMENTED DENSITY 

I Aril d201 lal) has shown that the Abel transformation of the 
augmented moment function of an anisotropic spherical sys- 
tem results in a similar integral transformation of the df as 
equation ( l2bl but with different powers on 7C and L 2 . This re- 
sult generalizes by means of the fractional calculus. The goal 
of this section is to establish them (see eqs.[36land[39ll for any 
pair of non-negative reals < fi < 

We start by considering to apply the integral operator of 
equation (0 to equation d2bl on ¥ or r 2 . In fact, we can es- 
tablish more general results. With 



d&dL 2( K s G(&,L 2 ) 



where the T* and r 2 dependencies of the integrable function 
G = G(£, L 2 ) are only through the two integrals of motion £ 
and L 2 (henceforth these trivial arguments of G will be sup- 
pressed for the sake of brevity), the Fubini theorem implies 



kfv^= ff dsdL^jyVw)] 

JJS>8 ,L 2 >0 



700(70 



+ ofA^)= ff d£dL 2 G o + f r2 V7O0(7O]. 

JJe>&o,L 2 >0 

T hrough direc t calculations that are basically identical to that 
of lAnl d201 lal appendix A) except for different arguments of 
the Euler integral for the beta function, we find that 



_ 0(7Q r* 
" r(A) J £+ za 

_ 0(7Q K s+A 
~ T(A) 2 A 
700(70 1 



dfiCF-fir 



L 2 

2{Q-S)- — 



B(A, s + 1), 



ofr 2 



r 2A+2 

r(A) 

0(70 r 



AR 



■(r 2 



r 2 Y 



(S-6) 

2A-2<jss+A 



HA) l 2a 

0(70 



j(2A+2 

B(A,s + 1), 



2Q¥ - £) - -7 



R 2 



T(A) 
0(70 



X 



R 2 ) A - 1 



dR 2 R 2s (r 2 

2U2-E) 

tfts+A 

B(A,s + 1). 



HA) 2 A ( x f-6) A 
Hence, we have established that 
Us 
2 A T(s 



i±iL_ ff dSdL 2 <K s+A G, 



,"K S+A G 



+£ At J s \ _ r 2A - 2 T(s + 1) ff 2 

Ob \—2)- r{s + A + l) JJ r d8dL L 2A • 

^K s+A Gd&dL 2 



+ c > (M r^ns + y rr 

oJ> V S')- 2 A T(s + A+l)JJ T 



(30a) 
(30b) 
(30c) 



which are valid for any s > - 1 and A > 0, provided that all 
integrals on the right-hand sides converge. 

We next find differentiations of the integral transform J^,, 
namely (here X = or r 2 ) 



s I 1 dfidirTO — G 

■L 2 sar t 2 



dX j 

The ^-factor for the .5 = case is due to 



(* >0) 



■Jt' 

Jo <3X 7f=o v 2r 2 7 



(31) 



8(70 = i8[T- — -fi] 

where o(jc) = 0'(x) is the Dirac delta. In addition, 
r2 [2r 2> T if£ = 



00 if S — —00 



Given that 



= 2; 



d<K_I? 
dr 2 ~ r 4 ' 



equation OTI ) suggests that for an integer n > and s > - 1 , 



PS. 



d£dL 2 70"G 



(n< s + 1) 



2 " <s, "/X' 

2T(1 + 5) dZ^G^ -,L 2 ) (n = s+l) 

{ Jo v 2H y 



(32a) 



WE 



d£dL 2 7C s "L 2 "G 



(« < 5 + 1) 



r(i + S ) 



rL 2 j2 

J dL 2 L 2 5+2 g(¥ - ^ ' i2 ) (» = * + !) 



(32b) 

Equations d30l l, d32l i and N = mo,o expressed as an integral 
transformation of the df as in equation d2"bl result in 



d" 



IX? -n) JJr 



(33) 



dBdL 2 -^ pr r(.B,L 2 ) (n<# 



2^£ 2 ^T(V-£ 2 ,L 2 ) (» = *)' 

ff dfidL^-"- 1 !, 2 "^,!, 2 ) (n<^) 
U? - n) JJ r 



XL 2 ,2 
dL^rfT-— ,L 2 ) (n = f) 

where n is again a non-negative integer and £ > |. Both 
equations further generalize from an integer n to a real fi < t; 
using fractional order derivatives, and it can also be shown 
that they are in fact valid for £ > if the extended definition 
in equation (0 is adopted. 
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In particular, to generalize equation (l33l l. we first find that 
+ r4 + r H( N \1 2 ^' 3 CCapat ^I ^F j 2 

[oh (— JJ = W(fn) JJ <&dL^rns,L 

T 

(35) 

for£ > i and A > 0, which follows equation (130b - The gener- 
alization of equation ( 1331 is arrived by applying equation d32i >. 
that is, for any reals < fi < £, and £ > 4 (the latter restriction 



that ^ > 5 will be dropped later in this section), 



(36) 



f 2 ^Vl^-3 rr 2 <k^ 1 

■JJr d8dL — ' 



r(^-yu) 



^(£,L 2 ) (/i<0 



2^1 r^' 



-3 



provided that the integrals converge. Equation 
now reduces to 



for £ = 



f 2 1+ ^ rr 

r(i-p)JJr 



dL 2 ^*?-— ,L 2 ) Oi=l) 



■ (37) 



Here setting \x = ^ —g results in equation d34t with n = given 

A— 1 1 

that i£, 5 N = + 6lp 5_f N. It is inferred that equation OH is 

in fact valid for £ > (n.b., < n < £, and so if < £, < \, 
then n = 0). 

A similar generalization of equation ( T34l i from an integer 
n to a real (cf., eq. |25T > and the extension of equation d36*l > 
to ^ > are possible although demonstrating them through 
direct calculations is comparatively nontrivial. Instead, we 
derive the generalization of equation (f34b followin g an indi- 
rect route. Let us first consider combining equation ( I30cl > with 
G - T, fi = s + 1 > and 1=1-6 where S = fx - |_//J, and 
equation ( f34b with n = and ^ = > 0, which results in 



; 7C W ^(£,L 2 ) 

r 



d&dL 



( 



r 4 — ) 



for /i > and < 6 < I. Next equation (I32bb indicates that 
for a non-negative integer n = |juj. However, 

(^d^) ^of,- 2 (^Sof/ 2N ) 



>J+4 | 



thanks to equation 



and consequently, we find that 

*L 2 r2 



(38a) 



This is also consistent with the case n = £ of equati on (l34b , 
again thanks to equation d25l ). That is to say, equation ( I38ab is 
actually valid for any fi > including integer values. 

Finally, consider a pplyi ng the integral operator in equation 
(01 on *P to equation ( I38al . as in 

where % > fi- The actual calculations is aided by an alternative 
expression for the right-hand side of equation ( I38ab 



(38b) 

It then follows that for < fi < £ 

(27T)2 

de(V-e)^ 1 dS(Q-STT[8,2r 2 (Q-S)] 
dS dQQP- Qf-"-\Q - SfT[&, 2r\Q - £)] 

^X^ff^K^^nS,^). (38c) 



(2r 2 



Equations ( I38ab and ( I38cb together, that is, 

mV'kfJM (39) 

j-^ Jf^dS dL 2 K^- l L 2 »T(&, L 2 ) (f > //) 

(? = /i) 



2K 



-4 dL 2 L 2 »T(y-±j,L 2 ) 



constitute the generalization of equation ( f34b from an integer 
n to a real which is valid for any pair of pi and £ with < 
fi < For < /i < £, < \, the indices transform (fi,t;) — > 
(| — £ ,| — jit) sends equation d39l to (T3~6b given equation (O. 
Equations (T3~6b and (|39l thus are both valid for any real pair fi 
and £ with < fi < 

In fact, both results and also equation[35]are different man- 
ifestations of the same result, that is to say, 



V2tt 



d & dL 2 ^ T ^ L2) ' 



2i-^ rr 
+ f + i)JJr U " 



(40) 



) 



which are valid for any real pair (A, £) such that /I + £ + | > 0. 
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4. MOMENT SEQUENCES AND AUGMENTED DENSITIES 

Consider the moment sequence of the df in (S,L 2 ) space 
restricted along ft = 0, given as in 



f df/ z/J^i/F; % r 2 ) (£ = 0, L 2 = 2r 2 ¥) 
Jo 



a) 



v Jo 



dF 7^(7; r 2 ) 



(£o = i 2 = °°) 



where 



#"(7; % r 2 ) = (27r)*TQ¥ - 7 2r 2 7). 



Then equations (f36b and (f39b indicate that 

^T^V^N) (p>\) 
l B ck^d^{r 2 ^) (Q<p<\) . 

U^o + f/(J) <£ = -M*q> 

In particular, if yu is a positive integer, this results in 



(41b) 



(42) 



N(Q,r 2 )dQ 



F„ = 







4~ — p f deCP - 0'-i(^)"[ r 2»N(e,r 2 )], 



where n = 1,2, That is to say, a set of fractional calculus 

chains of the AD directly determine the entire moment se- 
quences along a fixed sectional line in (£, L 2 ) space. In other 
words, the AD is similar to the moment generating function 
(or the characteristic function) for the df as a probability den- 
sity. With varying QP, r 2 ), the ft — lines eventually sweep 
the whole accessible (£, L 2 ) space, and thus N( V F, r 2 ) in prin- 
ciple uniquely determine the two-integral df, /(£, L 2 ). A few 
explicit inversion algorithms from N( V F, r 2 ) to f(S, L 2 ) are 
already available in the literature utilizing eith er the known 
inverse of named integral transforms (see e.g., iLynden-Belll 
Il962t lDejonghelll986t : IBaes & Van Hesd 120071) or comple x 
contour integrals (see e.g., Hunter & Qian 1993; Aril 1201 lah . 
Since the definition of the AD in equation (Q]i provides the 
explicit formula from f{&, L 2 ) to N( V F, r 2 ), the knowledge of 
N( V F, r 2 ) is therefore mathematically equivalent to knowing 
/(£, L 2 ). Once the potential T* = ?*(/-) is specified, the speci- 
fication of the AD thus completely determine a unique spher- 
ical dynamic system in equilibrium. Although this approach 
to the df /(£, L 2 ) through the AD N(*P, r 2 ) is advantageous 
as the observables constrain the AD more directly than the 
df, this procedure suffers a significant drawback in that the df 
recovered as such is indeed physical, that is, non-negative ev- 
erywhere in the all accessible subvolume of the phase-space 
- the "phase-space consistency", which is the subject of the 
reminder of this paper following the current chapter. 

Next, we consider what information on the physical prop- 
erties of the system is sufficient to specify a unique AD. First, 
we find from equation (l39l that the (augmented) velocity mo- 
ments of the even orders are related to the AD as in 



m Kn ?¥, r 2 ) = 



^ r 2„ + 2 \ r dr i) f efc ' (43) 



Here note that (^)t = T(k + ?)/ \[n. This is basically equa- 
ti on (13) oflDeionghe & M erritt ( 1992) - see a lso equation (8) 
of IBaes & Van Hesd ( 120071) equatio n (A2) of lVanHese etalj 
(2009), equation (5c) of lAnl(l2011bl) and so on. Equation 63) 
indicates that, given potential ^(r), specifying the AD com- 
pletely fixes every (in principle observable) non-vanishing ve- 
locity moment such that 

m Kn mr),r 2 } 



v 2k v 2n = 



N[¥(r),r 2 ] 

Conversely, equation d43l for (k, n) = (p + 1,0), that is, 

0+1. 



tn M+ 1 _o = 2^ + 1 ( j )+ j ^fy 1 * N at a fixed r reduces to 



MJr) = 



p\v, 



[T(r)F +1 f &q(f3t[qV{r)\r\ (So = 0) 
Jo 

o 



(44a) 



(So = -°°) 



where 



N[y(r)-(2,r 2 ] 
v(r) 



(44b) 



In other words, given the knowledges of the local density v(r) 
and the potential ^(r), the infinite set of the radial velocity 
moments in every order consists in the moment sequence of 
the AD considered as a distribution of ¥ - over the compact 
support if £ = or the half-open interval [0, oo) if £ = -oo 
- at fixed r. The problem is closely related to the HausdorfQ 
(for £o = 0) or the Stieltje^] (for So = -oo) moment prob- 
lems. With the infinite sequence of the radial velocity mo- 
ments as functions of r, the AD can then be uniquely deter- 
mined at least formally by such means as e.g., the HilberO 
basis or the Laplace and/or FourieiE3 transform (cf., the mo- 
ment generating function and the characteristic function) etc. 

The final information required for the full specification of 
the system is then the determination of the potential. Clearly 
the potential may be determined through the Poisson equation 
V 2 <E> = 4nGp, which under the spherical symmetry reduces to 



1 d/ 2 dT^__ 
r 2 dr\ 



dr / 



4/rGTv. 



(45) 



Hence if T = p(r)/v(r) is assumed to be constant, T(r) can 
be fixed by solving the ordinary differential equation on *F(r) 
that results from setting v = NOP, r 2 ) in equation d45l ). Alter- 
natively, from equation d43b . we deduce for k > 1 that 



dm k> „ 



d(r 2n+2 m Kn ) 
dr 2 



(2k- l)m k . 



(46a) 



(k- j)r 2n m k - Un+l . 



19 Felix Hausdorff (1868-1942) 

20 Thomas Joannes Stieltjes (1856-1894) 

21 David Hilbert (1862-1943) 

22 Jean Baptiste Joseph Fourier (1768-1830) 
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The total radial derivative of m^„ for k > 1 then results in 



dm, 



k.n 



2mt„ 



dr 



d\og(r 2n+2 m Kn ) 



d log r 2 



(n + 1) 



d*P dm k , n 
d7 d»F 



2(n + l)m tj „ - (2& - 



d¥ 

+ (2*- l)m H „- . (46b) 
dr 



With ¥ = *P(r) and m fcjB [T(r), r 2 ] = vufuf, this may be 
solved for d*P/dr if the required velocity moments as a func- 
tion of r are known. For the simplest case (k, ri) - (1,0), 
this reduces to the spherical (second-order steady-state) Jeans 
equation, 



1 d(vv 2 ) 2v 2 - 
v dr r 



dW 

dr' 



that is, the spherically-symmetric hydrostatic equilibrium 
equation with an anisotropic velocity dispersion tensor. 

5. NECESSARY CONDITION FOR SEPARABLE AUGMENTED 
DENSITIES 

In the following, we limit our concern to the cases for which 
the potential and the radius dependencies of the AD are mul- 
tiplicatively separable such that 



NOP, r 2 ) = P(^)R(r 2 ). 



(48) 



In addition to mathematical expediency, this assumption is 
also notable because under the separability assumption in 
equation d48l . the radius pa rt Rjr 2 ) of the AD alone can 
uniquely specify the so-called Binney anisotropy parameter, 



B(r) - I — = 1 — 

H 2v 2 r 2m lfi mr),r 2 ] 

1 d(r 2 m lfi ) 31ogmi i0 



(49) 



mifi dr 2 



8 log r 2 



'V(r),r 2 



such that dDeionghel Il986l lOian & Hunted fl99l 
Baes & Van Hesd 120071; lAnl 1201 lbt see also Ivan der Marel 
1994l as j? -1 being the integrating factor of the Jeans equation, 
i.e, eq.|47|> 



m = - 



dlog R(r 2 ) 



R(r 2 ) 

— ^ = exp 



>2f3(s) 



ds\. (50) 



dlogr 2 R(rQ 

Some application s are found e.g., in lBaes & Van Hesel (120071) 
while lAnl (1201 lbl) discusses further implications of the sepa- 
rability assumption. 

5.1. The radius part 

With a separable AD given by equation (|48T >. equation d39i > 
indicates that (hereafter x = r 2 ), 

^/ _i N) = ZffhQg) ■ + d/[x?R(x)] > (51) 
for ju < ^ whereas if 2 P > for £ > 5- Therefore, 

od/OfR) > (x > 0, fi > 0). (52) 
This is actually equivalent to the condition, 
d n [x n R(x)] 



R(n)(x) = 



dx" 



>0 (jc> 0,n = 0,1,2,...), (53) 



which is necessary for t he corr esponding df to be non- 
negative as noted by lAnl (1201 lbl) . It is clear that equation 
( 1521 implies equation (IB"3l as the latter is a restriction of the 
former fo r an in teger p — n. The opposite implication follows 
Corollary 12. 141 i: equation ( f53l > for a positive integer n implies 
equation d52l for p e [n — l,n] and thus equation d52l for 
p >0 follows equation d53l for all positive integers n. 

We find more equivalent statements of equation (T53t . First 
equation d25l > indicates that 



W = ^(4r[xR(x)] = (-l)V-^ 



where 



^(ur 1 ) 



(47) Hence equation ( 1531 is equivalent to 



(x 2 — )" [jcR(jc)] > (x > 0, 11 = 0, 1, 2, 



( -ir^> 

dw" 



(w>0, n = 0,1,2,...). 



(54) 
(55) 

(56) 
(57) 



Here the last is also equivalent to saying that the function l R(w) 
defined in equation ( l55l is a cm function of w. The Bern- 
stein theorem then indicates that 'Riw) is representable as the 
Laplace transformation of a non-negative function. That is to 
say, there exists a non-negative function <p{t) > of t > such 
that 'Riw) = Xf-»u»[0(O]- The inverse Laplace transformation 
4>{t) = Xzl+ffifw)] may be found using the Post-Widder for- 
mula ( fl91 l, which, thanks to equation d54l , reduces to 



hm —«(„)(-). 



Finally we find another equivalent necessary condition, 



lim 



1 d"[x"/?(x)] 



dx" 



x=//« 



> (f > 0). 



(58) 



(59) 



It is obvious that equation d53l implies equation ( |59l , pro- 
vided that it converges. The converse on the other hand fol- 
lows the Bernstein theorem and the Post-Widder formula. 
However, the conditional equivalence given the convergence 
of equation (l58l may also be inferred from equation ( 1271 ). 
By definition, equation d59l indicates that there exists a suffi- 
ciently large integer m such that R( n ){x) > for all v n > 3 m 
and x > 0. Corollary 12. 141 then suggests that R( m -\)(x) > 
for x > 0, and equation d53l l follows successive arguments on 
descending subscripts of R^ n) (x). 

5.2. The potential part 
IVan Hese et al.l (1201 lh proved that, given equation (J48j, 



p®0F)>o (jfc = o,...,L§-AJ) 



(60) 



where /3o is the limit of the anisotropy parameter at the center, 
is necessary for the df to be non-negative. We shall show that 
this generalizes incorporating fract iona l deri vatives . 

First, we generalize the result of lAnl d201 lal) to include ar- 
bitrary real order derivatives. This is trivial since the inverse 
Abel transform is just a particular fractional derivative as de- 
fined in equation (|4). If the AD is given as equation (|48| ). 
equation d36b reduces to 
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for n<^- Since R(x) > is trivially necessary, lj\o(x~ A R) > 
for x > and any A > u nles s R(x) = almost everywhere 
in x = r 2 e [0, °°) (Lemma [2.51 l. which will not be considered 
here. Consequently, equation ( f6Tb implies that 



°<5f/(5)< 



ia/p > o (ju < ^ + i). (62) 



With /I = 0, this indicates that ^<3/P > for any // < ~ - the 

condition for /i < is trivial because ^<Sy~' l P = g-fw/P while 
POP) > 0. For /I > on the other hand, equation ( I62i i implies 
that, if x~ A R(x) dx is integrable over x — 0, then id^P > 

for any /z < A + i and all accessible *P is necessary for the 
existence of a non-negative df. Alternatively, for a fixed y. > 
j, equation d62l suggests that ^d^P > is necessary for 

the df to be non-negative if there exists 3 A > /u - \ such that 

5f ( X ~ X R) i s well-defined. 

Equation d62b however is inconclusive whether tdy^P > 

is necessary for a non-negative df given R(x) ~ x~@ with 
P < 1 as x — » while this i s nece ssary if we were to extend 
the result of Van Hes e et al.l (1201 lb . For this, we first note that 
if /(f) is right-continuous at f = a, 



lim 



r h(t)6t 

1 £ r- 

+ Ja (t-a) l ~ 



lim hit) 



: /l(fl) (fl < fc). (63) 

This applied to the left-hand side of equation (f36t reduces to 



lim (| 



$of/ "(^-1/2) 



'7 



^-1/2; ^r(i-j/) 



(64a) 



where 77 < 1 and 

P„PF) = lim x"NCF,jc). 
Equation d36b overall then results in the formula 



'"^dO > 0, 



where 



lim L 2n T{&,L 2 ). 



(64b) 

(64c) 

(64d) 
(| - I?)" 

while maintainting < £ < | — 77. For 7/ = | - 77 on the other 
hand, the same limit is taken with fi = Therefore, this is 

valid for 77 < | - 77 and 77 < 1, provided that J Jr 2 (x^^N) is 
well-defined for £ < | - 77 (n.b., the integrability of the same 
for £ = I - 77 is actually 72of requ ired for its validity). Here the 
non-negativity of equation d64cb follows the n on-ne gativity of 
T{£>, L ). Of particular interests are equation d64ct for /u — 
and 1-77, 



For 77 < I — 77, this is derived with the limit £ 



p.m = 2l-%lr(i - iflSj^^W; 



2i"%ir(l -77) 



(65a) 
(65b) 



that is, explicit formulae for P,,(*F) and g^(*F) from each other. 
For a separable AD given as in equation d48l ). we have 



Therefore, equation (I64cb indicates that 

< R < 00 => ^<VP > (/j < § - 77). (67) 

That is to say, if there exists 3 77 < 1 such that R n is a (non- 
zero) positive finite constant, then td^P > for any V yu < 



|- 3 ?7. This actually encompasses equation d62b . which is seen 
as follows: If R^ is non-zero finite for 77 < 1, then we basically 
find that R ~ x-n as x -» 0. Hence & V*> converges for 

A < 1 - 77, and so if 7/ < A + 5, then 77 < | - 77. 
For example, with a constant anisotropy system given by 

R(x) = x~ ls ; Rp = l (J3<1) (68) 
the convergence condition reduces to 



y >(*) - J_ f C2L 



- s^" 1 ds _ F(l - j3 - A) 
~ xPY{\ -p) 



f A+f3 



< CO, 



which converges if < A < 1 -/3. It follows that equation 
indicates that + cVPCY) > 0for/j < < | -/J is necessary 



for the df to be non-negative whereas equation ( |64c| > suggests 
the same for p. < | — fi, 

6. SUFFICIENT CONDITIONS FOR PHASE-SPACE CONSISTENCY 
IN TERMS OF SEPARABLE AUGMENTED DENSITIES 



(66) 



Recently, IVan Hese et al.1 (120 121) derived the necessary and 
sufficient condition for the df with £ = to be non-negative, 
expressed in terms of the integro-differential constraints on 
the AD. They achieved this by reducing the problem to the 
Hausdorff moment problem, according to which the df is non- 
negative if and only if the moment sequence in equation ( 14 lab 
is a completely monotone sequencer^ Since the moment 
sequence are generated by the AD using equation d42l) . the 
monotone sequence condition is expressible in terms of finite 
differences of integro-differential operations on the AD. 

With a separable AD, they have derived a simpler sufficient 
(but not necessary) condition given as a union of conditions, 
each of which only involves the potential or the radius part 
separately but not together. Here we derive an alternative 
sufficient condition for a separable AD to be resulted from 
a non-negative df, wh ich turns out to be equivalent to that of 
IVan Hese et al.l d2012l) . The derivation here is based on the 
properties of cm functions and also uses the Laplace trans- 
form extensively. In this section, we only consider the case 
that £0 = and L 2 = 2r 2v F, that is, the df has a compact 
support and T(£ < 0, L 2 ) = 0. 

6.1. Inversion of a separable augmented density for the 
distribution function 

As it ha s been shown by dHunter & Oian 1993, see also 
lAnll201 la) , inverting equation d2bl for T{&, L 2 ) is formally 
equivalent to recovering the two-integral even df, T + (&, J 2 ) 
from the axisymmetric density v\*¥{R 2 ,z 2 ),R 2 ~\. The find- 
ings of the precedin g section together with the inversion of 
iLynden-Belll (119621) who utilized the Laplace transform for 
the latter problem suggest that the function (pit) defined by 
equation 081 must be directly related to the underlying df, 
T{&, L 2 ). We investigate this connection in the following. 

23 A sequence (00,01,02,...) is completely monotone if and only if 
(— 1)*A O; > for any non-negative integer pairs k and j. Here A is the 
forward finite difference operator defined such that A* +1 Qj = A*a J+ i - A fl; 
and A ci, = a,\ 
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Following Lynden-Bell ( 1962), we apply the Laplace trans- 
form on *F to equation I 



_£ [N( V F, r 2 )] = I d'Pe-^NCP 
Jo 



r 2 ) 



dfe- sT ^. (69a) 



£>0,L 2 >0 

The inner integral in the right-hand side reduces to 



I 

Jo 



Vf7<T 



, (69b) 



and consequently we find that 

V27T5 



X [N] 



Jo Jo 



d&e- s£ T(&,L 2 ). (69c) 



Substituting variables, t = \sl? and w = r 2 , this reduces to 



V 27T' w-* 



N(^, it;- 1 ) 



= £ 



r 

Jo 



dfie 



-.vfi. 



(70) 



If the AD is separable (eq. ED, then ur'NOF, uT 1 ) - 
P( y ¥)'R(w) where 'Riw) is as defined in equation (l55T > and so 
the left-hand side becomes 



r^O) 

(2*)s 



r<K f ) 

(2tt)2 



Here ^(s) is the Laplace transformation of 

!P(i)= X [PCP)] = I dTe-^ff) 



(71) 



(72) 



We have also used 'Riw) - £, t ^> w [<f>(f)]. Given that the in- 
verse Laplace transformation is unique, equating the right- 
hand sides of equations d70l i and (fTTb results in 



rdfie-^(fi,^). (73a) 

7T)I Jo V S J 



(2/r) 



Finally reinstating t = \ sL 2 leads to 

^(f)= Lineal ( 73b) 

The df is then recovered via the inverse Laplace transform, 

(74) 



(2n)i v 2 

6.2. Sufficient condition on a separable augmented density 

According to the Bernstein theorem, the df in equation d74i i 
is no n-negative if and only if the left-hand side of equation 
d73bt is a cm function of s > for all accessible values of L 2 . 
However Vis) defined in equation d72b is already cm since 
PQV) > 0. Hence Lemma IXTUl suggests that s?4>(sL 2 /2) is 
a cm function of s > for any L 2 > is in fact a sufficient 
condition for the non-negativity of the df. Equivalently, since 



d"[fi0(O] 



df 



t=sL 2 /2 1 2 I ds" [ ^ 2 I 



(75) 



the condition is also equivalent to that t?<f)(t) is cm. Unfor- 
tunately, this is too severe to be physically relevanQ which 
may be inferred from the constant anisotropy model given by 
equation d68l ). With this model, we find for/3 < 1 



d-y6)„ + 



0(f) 



1 



tPTQ. -p) 



where we have used 



lim 



i 



(76) 



(77) 



to find 0(f) using equation 

(13 -l)n 1 



r(i+z) 

The condition thus reduces to 



r(i -/?), 



>0 (f> 0,7i = 0,1,2,...), 



(78) 



for t > and all non-negative integers n, which cannot be 
satisfied for any constant /3 < 1 . 

Nevertheless, the preceding discussion extends to yield use- 
ful sufficient conditions. That is, for any fixed A, the condi- 
tions that 



„d*[^)] >() 
ds" 



(,y> 0,71 = 0,1,2,...), (79) 



(-1)"— — ^>0 (f> 0,n = 0,l,2, 
df" 



■) (80) 



are jointly sufficient to imply equation (I73bb being cm and 
consequently the non-negativity of the df. With increasing A, 
the constraint in equation ( 1791 tightens whereas the condition 
in equation ( f80b becomes strictly weaker. In other words, with 
a larger A, the smaller subset of functions PQ¥) will lead to 
s Ar P(s) being cm. At the same time if 0(f) satisfies equation 
(f80b for a fixed A = /lo, the same condition for any larger 
A > Aq automa tically holds. Both of these are easily inferred 
from Corollary 12.1 II 

6.2.1. the condition on R(x) equivalent to eq. \80$ 

Both conditions can also be translated into the direct con- 
straints on the behaviors of PQP) and R{r 2 ). For the radius 
part, we use that 0(f) may be given by equation d58l ). Note 
that the existence of 0(f) and the validity of equation d58l as 
well as its non-negativity, that is, 0(f) > for v f > are all 
necessary. Substituting equation d58l ) into the left-hand side 
of equation ( l80b results in 



( 1)B d"[fH0(f)] = ]jm (-1)" d" 



df" 



k\ df" 



lim 



(-1)" d n [x 3 2- A R {k) (x)] 



'(*) 



k\k> 



n+A- 



dx n 



x=t/k 



(81) 



Consequently, provided that this limit converges, equation 
( f80b is equivalent to insisting that there exists a sufficiently 
large integer m such that, for all integers v fc > 3 m, 



(~i)^{*i- 

^ ' dx n \ 



t d k [x k R(x)] 
dx k 



>0 (x> 0,« = 0,1,2,...). 

(82) 



24 If the Laplace transform of tp(t) exists, then 0(f) cannot diverges faster 
than t as t — > 0. Consequently, lim ( ^o t 3 ^ 2 </>(t) —> and thus ? 3 ' 2 0(f) cannot 
be cm because the limit suggests that ? 3 ' 2 0(r) should be negative or increasing 
in some interval t e (0, fo) where 3 to > 0. 
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That is to say, x^~ A R(^(x) being cm for all sufficiently large 
integers k is necessary and sufficient for <p(t) derived from the 
same R(x) to satisfy equation d80b . provided that the limit con- 
verges. In fact, equation ( f80b is equivalent to equation ( f82t 
for not only all sufficiently large inte gers b ut also all non- 
negative integers Ic, thanks to Theorem l2.15l which indicates 

that xi~ x R(m+i)(x) being cm implies x^~ A R( m )(x) being also 
cm. Succesive arguments with descending k then establish 
that x*~ x R(ja(x) being cm for all sufficiently large integers k 

implies that x^~ A R(^)(x) is a cm function for all non-negative 
integers Ic (the opposite implication is trivial). Note that the 
condition as stated in equation d82l> for all non-negativ e inte- 
gers k has already been noted by lVan Hese et al.l (120121) . 

6.2.2. the condition on -P( V P) equivalent to eq. \79\ 

The explicit constraints on P(?¥) resulted from equation 
d79l ) are derived by means of fractional calculus. We first 

find, from equations (O and ([171), that (n.b., ^Jp{0) = 
for £ > from CoroHary |2.71 i 



/Pis) = — 5 £ [PQ¥)1 = s» +1 £ Gf, ~°PQ¥)] 



(83) 



4*^5 J 



where n = \_A\ and 0<6-A-p.< 1. This indicates that 



> Q¥ > 0) 



^%0T= fa'PiO) = ■■■ = + (k A - l P(0) = 0, 



(84) 



(85) 



for A > is a sufficient condition for s Ar P(s) to be cm. 
Note, provided that P( V P) is right-continuous at *P = 0, that 

a_f<p ^(^) = ^ ^ c ^-' e 1- E3> which will thus be taken as 
granted. Consequently, equation d85l l for < A < 1 is essen- 
tially an empty condition. For A = 0, equation d84l i reduces to 
PQ¥) > 0. For a positive integer A = m + 1 on the other hand, 
the condition is equivalent to 



;>(»!+ 1) 



Q¥) > & P(0) = ■■■ = P m (Q) = 0. (86) 



For < 6 < 1, the condition in equation ( I85l l may also be 
replaced with the same bound ary co ndition as equation (l86l l. 
In particular, thanks to Lemma 127161 P (0) (0) = ■ • • = P w (0) = 
implies J<5p" +<5 / , (0) = for < 5 < 1. Consequently, it 
follows that for A > 1, 



p(0) (0) 



o(UJ-l) 



(0) 



(85a) 



actually implies equation (185t — if <5 = 0, they are identical. 
Therefore, equations d84l l and (I85at together also consist in a 
sufficient condition for s Ar P(s) t o be c m at a fixed A. The con- 
dition expressed with equation ( I85ab is useful because equa- 
tion ( fTTT i then indicates that equation ( f84l > is equivalent to 



1 d 1+ ^" f >f '^ ( " ) (0dg 

ep-e)« 



F(l -t5)dT 



> 



(87) 



where n is any non-negative integer not greater than A. 

Again, the joint condition of equations ( 184b and ( I85ab be- 
comes strictly stronger as A increases in accordance with the 
restriction on the complete monotonicity of s Ar P(s). This is 
seen using equation ( fT2b for < e < A under the condition of 



equation 



or ( I85al ). 



That is, J-V-PCF) > implies Jcl/POF) > for < £ < A. 
The similar implications of equation ( I85ab with descending /I 
are trivial. 

6.3. f/ze constant anisotropy model 

As an illustrative example, let us consider the constant 
anisotropy model with (3 < 1 (see Appendix 151 for the /3 — 1 
case) given by in equation ( f68b . Given equation d76l l, equa- 
tions ( f80l > and d82l now reduce to 



1)n d"[f 1-^(Q] = CS + A^jj" 1 > o 



df" 



r(l -j8) ^+"+-1-3/2 



(88a) 



,d J [it^, t) (i)] + {fi + A-\) + n 

(-1)" ■ „ ( ~ = (1 -ffi H ^Jn * 0. (88b) 



dx" 



JC /3+«+/i-3/2 



Thus /? + A > | and ft < 1 is sufficient for these to be satisfied. 

If A — m + 1 is a positive integer and /? > h-m, then equation 
d86l l i s sufficient for the existence of a non-negative df (cf., 
Ciotti & Morganti 2010). Our result furthermore implies for 
any r eal A > | that if | — A < /? < 1, then equations ( f84l > and 
HI constitutes a sufficient condition. 



With a fixed /3 < 1, this indicates that, if there exists 
3 A > \ - j3 such that equations (|84]i and (I85at hold, then the 
corresponding PQV) guarantees the phase-space consistency. 

. f d ^ n f n V c _< 3 - 



This also implies qcM^-P > for < f < A. while Sect 

indicates that, for the same system, qC^P > for v /i < | -jS 
is necessary for the df to be non-negative. It follows that, 

if N( v F,r 2 ) = r^PCF), then Jotp^P > is the necessary 
and sufficient condition for the non-negative df. In fact, here 
PQ¥) = pp(W) and T{&, L 2 ) = MS)Ir^ where Pp(?¥) and 
gp(S) are as defi ned i n equations ( I64bb and ( I64db with rj = ft, 
and so equation d65bt results in the inversion formula, 



+cfe^P(£) 



T(S,L 2 ) 



2§-%ir(l -/3)L 2 ^ 



NCP, r 2 ) = 



(89) 

This is sim ply the generaliz ed EddingtorF^l inversion formula 
(see e.g.. IE vans & Anll2006l) for constant anisotropy systems. 

That tb\t~Pp(?¥) > is necessary and sufficient for the ex- 
istence of a non-negative df is a trivial consequence of the 
inversion formula. 

7. FAMILY OF MONOTONIC ANISOTROPY PARAMETERS 

Consider the anisotropy par ameterized to be (cf, 
iBaes & Van Hesdl2007tlAnll201 lbl) 

Sir 2 ' 5 +B 2 r 2s 

m = 2, , (* > 0, r a > 0). (90a) 

If the spherical system is characterized by a separable AD as 
in equation ( l48l i. this follows the radial function (see eq. [30b 



R(x) = x~ pi (l + x s y^ where s£ = j3 2 -fa; 
<R(w) = uj- 1 R(w- 1 ) = ufi 2 -\l + w s y ( 

25 Sir Arthur Stanley Eddington (1882-1944) 



(90b) 
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where x = r 2 /rl (i.e., r. d = 1), which does not affect the 
following disc ussion. Note R(i)(x) > for x > restricts 
Pi , Pi ^ 1 • lAnl d201 lbh also provides a n elem entary proof that 
if < s < 1 and P\,Pi < 1, equation ( I90bb satisfies equation 
(T53T ). The same is deduced from the complete monotonicity 
ofR(w) for < s < 1 and p u p 2 < 1 (CorollaryEjJii, too. 

The situation for s > 1 however is inconclusive: on one 
hand, if p 2 — 1 > Pi, it is easy to show that ^'(w) < 
for < w s < (s — l)/(2 — Pi) and so the condition fails for 
s > 1 whereas lAnl (1201 lbl) on the other hand has found that 
the condition is met for all s > if — £ = (fi\ — Pz)l 's is non- 
negative integer. It appears that for a fixed s > 1, there exists 
a proper subset of parameter combinations P\,Pi < 1 that 
satisfies the necessary condition of equation 031 ). but we have 
not been able to establish the concrete criteria. Th e nec essary 
condition on the potential part discussed in Sect. 15.21 on the 
other hand is straightforward. That is, given R(x) of equation 
( 190bl ), the potential part P must satisfy ^dy^P > for any 
V /l < | - Pi in order for the df to be non-negative. Here also 
note Pi < 1 and thus ^ck A P > for any A< j. 

For < x < 1, the binomial expansion of equation ( I90bb 
and the subsequent term-by-term differentiation indicate that 



Zoo 
■ 
k=0 



k\ 



(91) 



k\ 

It follows equations ((58) and ( fTTT i that 



(-D*(a + f s ^' 



lim 



{\-p x +sk)l 



(92) 



k=o k\r(i -pi + sk) 



where E A pb (z) is the extended generalization of the Mittag- 

Lefflef 5 ! function E A nb i ntroduced by (Prabhakai1 ll971l see 

also lHaubold et al.ll201 ll) . Although the derivation here is es- 
sentially formal (see Appendix |A] for proper treatments) as 
in that we have not properly considered the issue of the con- 
vergence, the result is in fact valid given that Pi< 1 (for 
Pi = 1, see Appendix [Bb as is found in equation fllOOI i. Next 
we briefly detour to examine properties of generalized Mittag- 
Leffler functions necessary to derive sufficient conditions in 
Sect. [6] for the phase-spa ce co nsistency given that the radial 
part is given by equation J90bb . 

7.1. Generalized Mittag-Leffler function 

Definition 7.1 Let us consider a particular generalized hy- 
pergeometric function defined to be 



b (z) = £ 



Jp Un P k + b)k\ 



(P > 0). 



(93) 



Note that the Stirling^] approximation suggests 

r(«) .. i 



lim 

,i->oo Y(n + x) 



= lim 

n->oo (n + x) x 



(x > 0) 

1 (x = 0) 

oo (x < 0), 



26 Magnus Gustaf (Gosta) Mittag-Leffler (1846-1927) 

27 James Stirling (1692-1770) 



and so the ratio test for equation d93l with p > 
A + k zT{pk + b) 



lim 



k + 1 F(pk + b + p) 



= 



indicates that the infinite series for p > absolutely converges 
for all (finite) z- It follows that E A pb (z) with p > is an entire 
function of z- This is indeed a generalization of the Mittag- 
Leffler function since 

E l pb (z) = E lhh (z) ; 4, (z) = E p ,t(z) = E p (z) 

where E p (z) and E p j,(z) are t he classical Mi ttag-Leffler func- 
tion and its generalization by Wiman (1905). If p — 1 on the 
other hand, this reduces to the Kummec3 confluent huyper- 
geometric function of the first kind, that is, 



xF x (A-b-z) 



iFi(A-b-z) 
F(b) ■ 



Finally E pb (z) with A + is also a particular case of the 



Wright generalized hypergeometric function i W\ , 



1 



and also the Fo?0 H-function, 



(A, 1); 
(b,p); 



{1-/1,1} 
{0,1},{1 -b, P ) 



(94a) 



(94b) 



Next, the term-by-term integration indicates that for A > 



(95) 



where we have used Lemma 12.41 Together with the integral 
representation of the Mittag-Leffler function, 



E p , b (±z) 



~ f 



t''- b e'dt 
a f } +z ' 



and using (c.f., Lemma |2~4| and the binomial expansion) 

— T(k + A)(±z) k ti' A -PT(A) 
:z J= / — = 



<.=o 



k\tP k+ P 



(tp + zY 



this then leads to the integral representation 

f'^e'df 



=- r 



(96) 



Here the integral loop Q is the same as usual for the Mittag- 
Leffler function, that is, it starts and ends at '-oo', and loops 
around the circle \t\ — \z\ l ^ p in positive sense. This may also 
be independently proven using the Hankef^-loop integral 



1 1 f (0+) e^df 

F(i) ~ 2^i J_, ~ 



for the reciprocal gamma function (and also using the bino 
mial expansion), similarly to the classical case. Equation 



Ernst Eduard Kummer (1810-1893) 
Sir Edward Maitland Wright (1906-2005) 



29 

30 Charles Fox (1897-1977) 

31 Hermann Hankel (1839-1873) 
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implies the asymptotic expansion as z — > +°°, 



z A-b Q z 



P A W 



- z 



(-i)W 



1 



1 



(97) 



fe=0 



Jt!r(fe - pA- pk)z A+k T(b - pA) z 



If £ = -/I is a non-negative integer, the series in equation 
d93l terminates after the finite number of terms and thus re- 
duces to a polynomial on z - in particular, E° ph - I /T(b) is 
constant. In general, if £ = -A > 0, an alternative expression 
with the Fox H-function is also derived by separating the sum 
up to k = |_£J. That is, equation (l93l is alternatively given by 



Lf J / 

*=0 v 



f\ i-zf 



k) T(pk + b) 



U r (pk + b iap ) (k + m + iy. 



where b^ p = p(p + 1) + b and 6 = £ - |_£J- For < 6 < 1, the 
last infinite sum here results in 



r(i 



z* 



-.5) y 



= 2^2 



(i-<U),(i,i); 

(^,p),Gu + 2,l); 



-z 



{0,1} 



{6,1}, {0,1} 
-/i- 1, 1},{1 



The convergent integration path for the last Fox H-function 
with < 6 < 1 is always chosen such that it runs from c - ioo 
to c - ioo with < c < 1 - 6 whereas such straight paths 
do not exist for equation d94b| > with A < 0. Next, we find an 
extension of equation d95l > for a negative A = —g < 0, 



e ~ p %^ = z 



k=0 



(98) 



m-s) 



£ M (±z) - J] 



(±zr 



where = |_£J and <S = £ - /i. 
Finally we observe additional operational properties that 



d"^> fe (-z) 
dz" 



= (-l)W£^(-z), 



5f«X(-^) = Z 



(-l)m + 



fe=0 



1 



(l-^)„ + 



r(fo + pfe) r(n + k + i) 



«-] 



(99a) 
(99b) 



-j k\T{b -pn+ pk) 



for a non-negative integer n. The last holds given that (1 - 
A)* * 0. In addtion, using (A)\{A + k) = (A) + k+1 = A(A + 1)+, 
we also find that 



d[z^(-z)] 
dz 



k=0 



k+A-l 



k\T(b + pk) 



(99c) 



7.2. Sufficient conditions for the phase-space consistency of 
eg. MM with < s < 1 

Now we consider sufficient conditions on the AD to guaran- 
tee the phase-space cons isten cy (Sect. [6]l with the radial func- 
tion giv en b y equation J90bl > with < s < 1 and £q = 0. 
In Sect. 16.31 we have argued that for P\ = A < 1, if there 
exists 3 A > | - Pi = \ - p 2 such that > and 

P(0) = • ■ ■ = p(UJ-D = o, then the df with fi = recovered 
from the particular PQ¥) and R(x) is non-negative everywhere. 
This follows from the fact that f^'V(0 = ^ -,HJ /T(l-/?)iscm 
if and only if A > |-j8. As with t ^tpff) = fi j (-f) 
for j8i < 1, the discussion in Sect. 16.31 on sufficient condi- 
tions for constant ft separable AD can carry over here essen- 
tially verbatim if we can establish the set of results regarding 
the complete monotonicity of the generalized Mittag-Leffler 
functions of the form of f'E ( s Y _„ (-f). 

We first note that the leading term of E ( , „ (-t s ) for t ~ 
is given by the positive constant 



l 



r(i -A) 



> o 



for A < 1, which indicates that fE l sl _^{-f) for A < 1, 
5 > and a > must be increasing in some interval (0, c) 
where 3 c > 0. On the other hand, equation d97b suggests that 



lim i 



-Pi fit l-t p> i 



1 



m -j3 2 ) 



> 



we have fE ( , . (-t p ) 



are positive finite for /?2 < 1. It follows that, as f — > +oo, 

+oo for fii,/3 2 < I, s > 0, and 

a > Pi That is to say, if a > p 2 - Pu then fE l sX A-t p ) 
with for A ,/?2 < 1 and 5 > must be increasing in some non- 
empty subintervals of (c, oo) where 3 c > 0. Together, these 
observations imply that fE^ t (-t*) for Pi < 1, p 2 < 1, 

and s > cannot be cm if a > min(0,jS2 - y£?i ). Although it 
is tempting to hypothesize by analogy to the constant p case 

such that fE l sl A-f) for p x < l,fi 2 < 1, andO < 5 < 1 are 
cm for a < min(0,jS2 - Pi), we have only been able to prove 
this under a restriction that p 2 > Pi or p 2 < 1 - s while for 
1 - s < p 2 < P\ < 1 we only manage to find a more restrictive 
condition a < -s (n.b., -s < -s + 1 - P\ < p 2 - P\ < 0) for 

r 

the complete monotonicity of fE s yai—f). 

We next apply the similar discussion as in Sect. 16.31 and 
find that for R(x) given by equation ( I90bb with fli < l,p 2 < 1, 
and < s < 1, if the re ex ists 3 A > I - min(y6i,yS2) such 
that equations d84l i and (I85al > hold for P(90, then the resulting 
AD, Pi^Rir 2 ) guarantees the existence of a non-negative df, 
unless \ - s < p 2 < P\ < \. For (s,p 2 ) - (1, 1), which results 
in E'-f^i-t) = !Fi(l -A; 1 - A;-0 = e-'^iCO; 1 -Put) = 
exp(-f)/r(l — Pi), the condition for a positive integer A — m + 
1 > \ ~ Pi (note min [ ffi,/?2l = Pi < 1 = A) repr oduces that 
of ICTotti & Moreanfil dioioi) . for the generalized ICuddefordl 
system to result from a non-negative df. If 1 - s < p 2 < 
Pi < 1 on the other hand, we at this point only find a slightly 
restrictive sufficient condition with 3 A > \—(fii—s) > \~p 2 > 

| -A > 3 ( n - b -<A -* < 1 -s <A <A < !)■ 
In the rest of this section, we proceed to prove that 

fE e sl _ pi (-f) where a = min(0,A - A) and ^ = A - A for 
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y8i < 1, j8 2 < 1, and < s < 1 (but not 1 - s < /3 2 < (3i < 1) 
is cm. First, if /Ji = #> < 1, then ^^(-f 1 ) = 1/T(1 - jSi) 
and so is trivial (see Sect. 16.3b . Next note, if we can prove that 
^mn(0A) E j ft (_ z ) i s cm for b > 0, b > pA, and < p < 1, then 

the desired result follows Corollary |2.11l In the following, we 
prove the complete monotonicity of E pb {-z) for < pA < b 

and z A E pb (-z) for b > and A < 0. The further restriction, 
/?2 < 1 — i (i.e., + p < b) on the latter case meanwhile 
occurs naturally. We first introduce a lemma, 

Lemma 1.2 If < p < 1, b > 0, and b > pA, then E pb (-z) > 
is non-negative for all z > 0. 

This formalizes the fact that fi(w) for < s < 1 and fii ,02 < 1 
is a cm function of w > 0. In general, for b > 0, 



dte~ wt t pk+b - 1 



/>=() 



k\wP k+b w h 



(100) 



By Corollary |2.12l this is a cm function of u; > forO < p < 1 
either if £ > and ^ < or if b - pA > and A > 0. Then 
from the Bernstein theorem, if0<p<l,fe>0, and b > pA, 
then t b ^E A b (-t p ) > for t > and so E A b (-z) > for z > 0. 
The first half of the desired result is now trivial, that is, 

Theorem 7.3 IfO < p < 1 ami < pA < b, then E pb (-z) is a 
cm function of z > 0. 

This follows equation (I99ab. Note that if b > pA, then b+pn> 
A + pn and thus Lemma I7T21 to gether wit h (A )* > for A > 
completes the proof. As noted, Theorem l7 . 3 1 implies 



Corollary 7.4 For &q — and R(x) given by equation A90b[ 
with < s < 1 and Pi < fii < 1, if there exists 3 A > j - f>i 
such that > and P(0) = ■■■ = p(W-D (0 ) = o, then 

the df inverted from P( v F)P(r 2 ) is non-negative. 
This actually extends to fii < 02 < 1 (Sect. I6.3l l. Also note 

that if P(0) = • • ■ = p(li-^)(0) = 0, then Jo^'P > is the 
necessary and sufficient condition for the phase-space consis- 
tency given &o — and R(x) with < s < 1 and fJi < f3i < 1. 
For the second half, we first find 

Theorem 7.5 If < p < I, b > 0, and £ > 0, f/ien 
z) is a cm function of z > 0. 

Corollary 7.6 For £q — ant/ R(x) in equation H90b[ with 
< s < 1 and fa < ySi < 1, if there exists 3 A> | — jSj + sn 
where n — \(J3i —/32)/s~\ such that qO\i A P > one/ P(0) = 
• • ■ = j p(W-i)( ) = 0, f/zen f/ze df inverted from P( v F)/?(r 2 ) is 
non-negafive. 

If £ = // is a non-negative integer, this is trivial since 

7 -(M-k) 



(101) 



with every coefficient being positive. Next, equation (I99al > for 
£ = — f < and n - \tf \ results in 

d^E'M-z) 



dz^ 



= (1 - e) 



i-z), 



(102a) 



where < e = — t; < 1. Now equation ( |99b| i indicates that 

(102b) 

which is consistent with equation ( fT2l . If £ > 0, this reduces 
to (note then that \f \ > 1) 



(1 - e)+ r i 

' j d«(l- «) m - 1 ^ fc4 . Bm (-w). (102c) 



(tel " 1) 



j p,i+prfi v 



If e = (i.e., r^l = f), then (1 - e)^ = f ! and & pJb+pg = 
1 /F(/7 + p%) > 0, and so this is just equation (1 1 1 b . In general, 
this implies that z~^E~f b (-z) with < p < 1, b > 0, and 
£ > is cm since 



n r\ 



X 



, dw (1 - ufHzu) = I da(l -u) k u"f"\zu), (102d) 

di" 

while Theorem 17.31 indicates that f(z) = E e pb+p ^(-z) is cm 
given b + p\%~\ — pe = b + pi; > 0. Finally, we are able to prove 

Theorem 7.7 IfO < p < 1, £ > 0, b > 0, and b > p{\ - 

then z~^ E pb {—z) is a cm function of s > 0. 

Corollary 7.8 For Sd = and R(x) given by equation A90bi 
with < s < 1, p2 < 01 < 1, and 02 < l—S, if there exists 3 A > 
such that *dv A P > and P(0) = • • • = ^"''(O) = 0, 
f/ien the df inverted from P(*F)P(r 2 ) is non-negative. 

If £ is a positive integer, this is t he sa me as Theorem 17. 5 1 For 
general cases, we note equation d99c| ) results in 



d[ Z -%%-Z)\ ZE'Ji-z) _ tZ-^E-f\-z) 



dz 



T 2-e 



(103) 

where [£-11 = ^1 - 1 and < e = \%\ - % < 1. Theo- 
rem [73] indicates that if0<p<l,fc>0, and ^ > 1, then 
z U-£J g~f (_ z ) anc [ subsequently z~^ +v *E ~£(—z) are cm. The- 
orem l7.3l on the other hand suggests that if < p < 1, £ < 1, 
and b > p(l — ^), then E pb {—z) is cm. Hence, if < p < 1, 

£ > 0, b > 0, and b > p{\ - g), the derivative of z~ ( E~ p \{-z) 
is given by a cm function multiplied by a negative constant. It 
follows Lemma |2~T0l 2 that z^E p {(-z) for < p < 1, ^ > 0, 
> and > p(l - is a cm function of z > 0. 

8. SUMMARY 

We have shown that the fractional calculus operations 
(eqs.[3]H] and [5]) applied to the bivariate augmented density 
(eq. []]) result in a set of the integral transformations of the 
two-integral distribution function (eqs. [36l[39l andl40li. Equa- 
tion (1401 with A + g + j =0 indicates that the set of frac- 
tional calculus operations on the augmented density N( V F, r 2 ) 
listed in equation d42l provides with the complete moment se- 
quence of the distribution function along "7C(£, L 2 ; l F, r 2 ) = 
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as shown in equation ( 141 ab . We infer from this that the aug- 
mented density that ensures the non-negativity of the distri- 
bution function may be deduced by analogy to the classica l 
moment problem in probability theory (IVan Hese et al.l20l2h . 
We have also found that equation ( f4Ub for a non-negative in- 
teger A > and £ = consists in the complete moment se- 
quence of the augmented density at a fixed r considered as 
a probability density on *F - which is possible because the 
augmented density is also non-negative in all accessible r and 
*P. Comparing this sequence to the velocity moments result- 
ing from the given distribution function (eq. |2bl , we deduce 
that the augmented density (and subsequently the distribution 
function) is uniquely specified given the potential ^(r) and the 
density profile v(r) once the infinite set of the radial velocity 
moments in every order (equivalently the complete radial ve- 
loci ty distribution) as a functi on of the radius are available 
(cf. JDeionghe & Memtll992l) . 

Given A + £ + | > 0, all the integrands in the right-hand 
sides of equation (l40l > are non-negative because the distri- 
bution function T(&, L 2 ) must be non-negative in the whole 
accessible subspace volume T. This non-negativity implies 
that it is necessary for the integro-differential operations on 
the augmented density N(*F, r 2 ) given in the left-hand side of 
equation d40b to be also non-negative, provided that the in- 
tegrals involved in their definitions are all convergent. This 
introduces the set of necessary conditions on the augmented 
density for the non-negativity of the distribution function. If 
the augmented density is multiplicatively separable into func- 
tions of the potential and the radius dependencies li ke equa - 
tion d48l ), this results in the condition stated by d201 lbl) . 
that is, equation (l53l for the radius part of the augmented den- 
sity. We have also discovered a few equivalent statements of 
this condition, notably equation d59l l and the function r R(w) 
defined in equation 05l ) being completely monotonic and so 
on. The same argument for the potential part of a separable 
augmented density on the other h and recovers the conditions 
derived by IVan Hese etai] (120111) and |Xn| (1201 Id) . They are 
further generalized with fractional calculus to indicate that: 
g^oly^P > for all accessible is necessary if p < i or there 

exists 3 A > yu - i such that ^j- 1 \r~ 2/i R(r 2 yi is well-defined or 

3 jS < | - fi such that lim,.2^ + r 2 ^R(r 2 ) is non-zero and finite. 

With separable augmented densities, the distribution func- 
tion may be inverted from the augmented density by means of 
the inverse Laplace transform as in equation ( T74b . The non- 
negativity of the distribution function corresponding to a sep- 
arable augmented density is guaranteed if the Laplace trans- 
formation of the distribution function given in equation ( 173b| > 
is a complete monotonic function of s > for any L 2 > 0. 
We have shown from this that the set of joint conditions com- 
posed of equation ( l82b with all non-negative integer pairs n 
and k for the radius part R (r 2 ) of the augemented density and 
equations (l84b and J85at for the potential part PQP) of the 
same is sufficient to imply the non-negativity of the corre- 
sponding distribution function. This last set of suffici ent con- 
ditions is equivalent to that of I Van Hese et al.l (120121) . which 
was derived from the argument following the application of 
the Hausdorff moment problem. 

This manuscript is basically an extended version of lAn et all d2012h . 



APPENDIX 

A. PROPER DERIVATIONS OF <j>(f) EST Eq. {92) 

Let us define with R(x) given by equation d90b| > so that 



R 



1 + X s 



T„ EE (1 + X S )"a n , 



y = x 



(Ala) 
(Alb) 



Using equation d26j >, we find the recursion formula for a„{u), 
1 d(x n+1 Ra„) dlog(x" +1 fl) dlogwda,, 

Gn + l - — T7 : = — d-n + U- 



x"R 



Ax 



d log x 



d\ogx du 



= [n + 1 - B\ + (B\ - B 2 )u] a„ + su(\ - «)— A (A2a) 

du 



Given that 
da„ dy d 
du du dy 



(1 + y)" 



l 



dr„ 



(l+y)"- 2 dy (l+y)"- 1 ' 

(A2b) 

the recursion formula for r n (y), 



t„+i = (1 +y)" +1 a n+ i 



n + 1 - Pi + 



l+y 



(1 +y)" +1 a„ + sy(l + y) 



(A2c) 
n-l^n 

du 



dr„ 



= [n + 1 -Pi + (n + 1 -p 2 - pn)y]r n + sy(l + y)—, 

dy 

also follows. Both recursion formulae imply that a„{u) and 
T n (y) are an (at most) n-th order polynomial of their respec- 
tive arguments, u and y (note oq — tq — 1 by definition). 
Subsequently, if we let 



then 



n (l n k x 



A-Pi 



=« (l +x s ¥ +k (i+x s y + 



(A3a) 



(A3b) 



In addition, given that t„ = (1 +y)"a„, the binomial expansion 
and the subsequent rearrangement of the double sum 



t„ = (1 + y)"a„ = a n , k y k (l + yf k 



k=0 



n m 



n n-k 

EE 

k=Q 7=0 



n — k\ 



k+j 



(n-k)\ 



a n,ky 



a n,uy 



(A4a) 



leads to the relation between the two sets of coefficients, 

1 v-i»i (n-k)l . .. ... 

tn,m = 7 T7 > .. n ; ~ a n,k- (A4b) 

(n - my. *—v<=q f 



(n - m)\ ^—'k=o (m - k)\ ' 
We note that if we define the associated coefficient sets, 



tfi,k — I , I tn,k s 



a n ,k = (-1) W\ a r i 



(A4c) 



the relation reduces to the standard binomial transform, 

k /, \ k 



tn,k 



m=0 



m 



a n ,k 



m _0 \ / 



which is known to be involutionary. 



(A4d) 
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The expression for the coefficie nts is found using the bino- 
mial series expansion of equation ( IA3bb for < x < 1, 



= 2j a ">* 2j •) X 

k=0 j=0 



CO m 



m=0 A-0 



(A5a) 



(m-k)l 



where we have used 5„ t = if k > n. All the Pochham- 
mer symbols without any directional specification hereafter 
are interpreted to represent the rising product, i.e., (a) n — (a)^. 
Match ing the coefficients for the same power of x in equations 
dA5ab and (EQi leads to [n.b., (£) m (£ + m\_ m = (Ok] 



(l~p 1 +sk) n =Y 

X—J»(=0 



k (-iy n k\ a, v 



(k-m)\(On 



(A5b) 



Although its derivation assumed < x < 1, this is valid re- 
gardless. The right-hand side is in the form of the binomial 
transformation and thus by its involutionary inversion 



(-1)* 



*=o kl(m - k) 



T (l -ySi +Sk) n 



(A5c) 



That is, is the £-th order forward finite difference of (1 - 
A +sx)„ = JTJ_j (J — A + sx) at x = 0. Since (1 -A + sx) n is an 

«-th order polynomial of x, we have A*(l - A + sx) n \ x= p = 
if k > n. The formula for t„ is found from equation dA5cl ), 



(-D" 



-V. . 
Z — im={) 

- (£)n+l > . . — 7 
x — Jm=0 mj! 



(o;(«+^- m (i-A + *< 



m!(fe-m)! v ' /mv ' ' ' /<; -" 1 
(-1)'" 



(A6a) 



"=0 m\(k - m)l (( + m)i +n - k 
using equation (IA4bt and the Chtf^l-Vandermondd^l identity 
Z" k= o{"k)(sMt)n-k = (s + t) n , (A6b) 
or equivalently the Gaus£3 hypergeometric identity 
J" ( -ifl n )Vk = Y n (-nh(b) k = (c-b) n 

^k=Q V \k) ( C )k *~>k=0 k\(c) k (C) n 

For the 5=1 case, from equations (lA5cl >. dA6cl >. and 
(1 - AM 1 ~Pi +k)„ = (1 ~P\ +n)k (A7a) 

we find that (A < 1) 



(l-A)n m!(l-A). 



(A -/3i)„,(-n) m _ /_^j«l n l 0^2 ~A)» 



m/ (1 - P\)„ 



tn.k 



(i-A)« 



£(-D" 



m=0 



fc\ (A~A) m _ (1 ~A)t 
m) (1 -j6i)„, (1 -A)* 



(A7b) 



(A7c) 



32 Zhu Shijie (1270-1330) 

33 Alexandre-Theophile Vandermonde (1735-1796) 

34 Johann Carl Friedrich Gauss (1777-1855) 



Here £ = A — A since s = 1. This is notable as it indicates 
that t„, k > 0, and with (1 - P\) k (l - P\ + = (1 - Pi)„ that 

r n (y) = Z(3 (1 -A + &)«-*(! -Pi)ky k > (A7d) 



R {n )(x) = 



(1 +xf2-Pi+n 



Pi+k) n - k (l -p 2 )kX k ~ Pl 
(A7e) 



for any non-negative integer n and all x — y > 0. 
Next, we consider the Mellin^ transform for < z < A 



tp(z) - M 



I 

Jo 



du(l — u) 



-f 

Jo 



<>'n 



(1 + yY 



(A8a) 



-I 



r( z + fe)T(j - Z ) . 



r( Z )r(A-z) 
r(i) 



z 

A=0 



A v (l — B\ + sx)„\ 



This simplifies for A = C, utilizing the Newtor0 series 

001 



(A8b) 



If /(jc) is an n-th order polynomial, the formula is exact after 
the summation up to k — n. Since A*(l + sx)„\ x= o = for 
k > n with the n-th order polynomial (1 — /3\ + sx) n , 

(1 - p 1 - sz) n = 2I =0 ( "^ (Z) * A*(1 - A + «)»U- (A8c) 
and therefore with A = £ > z > 0, 

r(o^(z) = r(z)r(f - z )(i - a - «)„ . (A8d) 

By means of the inverse Mellin transformation, /?(„) = /Ja,, is 
then expressible to be a Mellin-BarnesFl type integral 



R(n)(x) 



i r dz 

2m x@ l J c x s: 



dzn Z )T^-z)(l-pl-Sz)n 



no 



(A8e) 



Although this is actually reducible to algebraic functions on 
x 1 as in equation (|A3bJ, it is also the iFoxl H-function (H^ 
in particular) and further reduces to the MeiieiPl G-function 
(G '"i j) and the hypergeometric function („+iF„), the last of 
which would be formally equivalent to equation d9lT l. 
The function 0(f) is found from equation I 



1 Cdz 

~ 2mt^J c t sz 

- 1 r^£ 



roone-z) lim (i-A 



SZ\ W 6 



no 
nz)n(-z) 



t^no 



c r~noni-Pi-sz) 

{0,1},^!, 5 



(A9a) 



where we have used equation ( |77] i. IfO < ( < (1 - Pi)/s 
(n.b., ^ = (A - P\)l s and A ^ 1)> the convergent integration 

35 Robert Hjalmar Mellin (1854-1933) 

36 Sir Isaac Newton (1642-1727) 

37 Ernest William Barnes (1874-1953) 

38 Cornells Simon Meijer (1904-1974) 
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path C may be chosen such that z - c — z'oo to z - c + z'oo 
with < c < Provided that £ is neither zero nor a negative 
integer, this is still valid but the integration path should rather 
be chosen to separate the poles of T(z) from those of F(£ - z). 

In fact, (pit) may alternatively be found in terms of an infi- 
nite serie s for a ny (, In partic ular, eq uation d58l after inserting 
equation (lA5cb into equation dA3bb results in 



n k 



R(n) - 



yy i-mO* v 



sk-fii 



(i +x s Y +k 



(A9b) 



n k 



0(0 



Hm ZZ ( ~ lfl ° k 11 ~ 0: + sq) " 



k=0 q=0 



(1 +X s Y +k 

(1 -ySi + sq) n 



x=t/n 



Zoo 



q\(k - q)\ n^oo n \n sk -^(l + t s /n s Y +k 

(-D k (Ok 



t sk- 



=o k\T{\ -pi + sk) 



1 



(f, i); 
(l - A, j); 



-r 



where we have used equation ( ITTl i and 



lim — 

n— >oo yi a 



1 (a = 0) 
(a > 0) 



(A9c) 



(A9d) 



For 5 = 1, we have 0(f) = r 13 ' i A(#> 1 ~P\;-t) where 
\F\{a;b;x) = iFi(a;b; x)/T(b) is the regularized hypergeo- 
metric function. The non-negativity of <p(t) > for t > 
and s = 1 is explicitly shown by the Kummer hypergeometric 
transformation resulting in (f>(t) = r^'e~'iFi(l — /fo 1 — j3\\f). 
That is to say, from equation d92l in general 



exp(f 5 )<X0= J] 



j.s{k+m) 



, -J jfc!m!r(l -fli + s/fc) 



n=0 *=0 W 



(0* 



r(i-j8i + sfe) 

(A10a) 

The inner sum for s — 1 simplifies, from equation (IA6cl ), to 



H3 ! 



PCI -jSi +Jfc) rfl-jSi+n)' 
and therefore for f > 

(1 ~Pz)n 



e 



■f>0, 



(AlOb) 



(AlOc) 



which is non-negative for all pairs (fi\,p<£) given that /8i,/?2 ^ 
1 as every coefficient of the series is then non-negative as well. 

B. THEySi = 1 CASES 
B.l. A proof of eq. 
First, we note a trivial result, 

Lemma B.l For c, A > 0, 



A f x 11 Ax 
Jo 

Next, it follows that 



lim c A 



1. 
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Theorem B.2 for F(x) = f(x) - I where I = lim A -^o+ fix), 

lim A f jC* -1 |F(*)| dx = (c> 0) 
<^<> + Jo 

proof. First by the definition o/lirn t _ >0 + f(x), we find that for 
any v e > 0, there exists 3 5 > such that, if < x < 5, then 
\f(x) - £\ = \F(x)\ < e. Now ifO < c < 6, then for any A>0 

0< f x^lFix^dx < e f x^'dx. 
Jo Jo 

/f c > 5 > on the other hand, 

fC ffi r*c 

0< I z* -1 |F(jc)| dx = I x^'lFWIdx-H I z* _1 |F(*)|dJt 
Jo Jo Js 



< e 



f x^dx + sup[|F(jc)|] f x^'dx. 

Jo (S,c) JS 



Note here that x A 1 dx is finite. Consequently, provided that 
f(x) is bounded in (0, c), we find from both cases that 



< lim A f x^'|FW|dx 
^->0+ Jo 



< e 

where c > ant/ we have used Lemma \B. 1\ q.e.d. 
It immediately follows that 

Corollary B.3 for F(x) = f(x) - I where I = lim A ^o+ fix), 
lim A [ x^'F(x) dx = (c> 0), 

A-* * Jo 



lim A f 

-i-o + Jo 



x A ~ 1 fix)dx = £ (c>0). 



Equation (|63T > trivially follows this with the change of integra- 
tion variable x-t-a. Formally this is interpreted to be 



lim Ax 1 ' 1 = 6(x) ; lim l - 

A^0* a ^\- X C T(1 - a) 



= 6(x) (Bl) 



where 5(x) is the DiracPl delta, provided that /(x) is right- 
continuous. 

B.2. The p — 1 constant anisotropy model 
Let us consider the df given by 



, fiS) 6(L 2 ) 



(B2a) 



where f(&) is an arbitrary function of £. This df corresponds 
to the spherical system entirely built by radial orbits, that is, 
the /3 = 1 constant anisotropy model. Given that 7C(L 2 = 0) = 
2( V F - S), the corresponding AD is found to be 



N( l F, r) 



V2( v F - 8) 

which is separable as in equation (l48b with 

39 Paul Adrien Maurice Dirac (1902-1984) 



(B2c) 
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The AD is inverted to the df using the fractional derivative, 



/(£) = ^-P{&) > 0, 



(B2d) 



whose non-negativity is also the necessary and sufficient con- 
dition for the phase-space cons isten cy. Note that this is con- 
sistent with the results of Sect. 16.31 applicable for /3 < 1 as is 
R(x) here the natural limit of the constant anisotropy model in 
equation d68b to ft = 1 . 
Furthermore, we find for A - n + 5 > and n - \_A\ that 



(x-y) A l dy 



V 

x ,,A-\ 



,1+1 



y Ay _ x"F(A) 
{x-yf ~ n\ ' 



(B3b) 



i-6 A _i _ T(A) d n+l x n 



n\ dx n+1 



0, (B3c) 



while ^j- °x 1 = q3 a °x 1 = x l . Hence, R(x) = x 1 satisfies 
the necessary condition in equation d53l ). Moreover, equations 
(HH), ([39>, and (@0j still hold with non-trivial cases indicating 



So 



(B4) 



whose non-negativity for V ju < \ is t he same necessary condi- 
tion for PQ¥) discussed in Sect. 15.21 

From R(x) = x , we find that 'Riw) = 1 and its in- 
verse Laplace transformation at least formally is given by 
0(f) = 5(f). Although equation ( f80b strictly is then trivial as 
5(f) = for f > 0, this interpretation of equation < T80b seems 
improper considering that the Dirac delta is not differentiable 

at t = 0. Equation (T82t on the other hand reduces to x^~ A 
being cm since R(Q)(x) = R{x) = x~ l and R(„)(x) = for any 
positive integer n. The sufficient condition following this, that 
is, equations (l84l and d85al i for 3 A > \ is in fact a proper one, 
as is the natural limiting case of the constant anisotropy model 



for p = 1 . It appears that for R ~ x as x ~ (and lim„ 



being nonzero finite), we may consider (f>(f) ~ t as f ~ for 
the purpose of applying equation (f80b . 



B.3. Equation A90bi with fix = 1 

The discussion in Sect. on necessa ry co nditions is valid 
inclusively for f3\ < 1. That is. equation ( I90bb with ffi = 1 still 



requires to satisfy equation (1531 - if < p < 1, this is auto- 
matically met - in order for the df to be non-negative whereas 

i 

the potential dependent part is restricted to be ^oly J P > for 
the phase-space consistency. 

The complication arises however for /3\ — 1 in regards 
to sufficient conditions discussed in the preceding section. 
The main difficulty is due to the fact that lim^o^M = 
lim^oo 'Riw) = 1 is non-zero. Whilst this would indicate 
(f> ~ r l for f ~ 0, the particular behavior is incompatible 
with the convergence of the integral. The formal solution fol- 
lows adopting equation (IB lb . In addition, the limit of equa- 
tion d59l with R = x is identically zero for any x > and so 
the function <p(t) defin ed vi a the formal limit of equation d5"8l l 
with R(x) in equation J90bb takes the same value as that with 
"R(x) - x" 1 " for all f > (that is to say, the Post-Widder 
formula is technically valid). In other words, the function 
(f>(t) derived in equation [92] with fi\ — 1 is in fact the inverse 
Laplace transform of "^(w) - 1" and the 'true' inverse trans- 
formation of Hiw) with /3\ = 1 is given by <p(t) + 6(t). For 



example, since 1 /r(0) = 0, the k — term in the power series 
defining the generalized Mittag-Leffler function E Q does not 
contribute. Hence, equation ( 11 001 ) can in fact be well-defined 
for the b — case too. In particular, 



d>e-«r'£j (-O 



Z 



*=i k\wP k 



i + -L, 

wP 



L. 



Since (1 + w^ p Y ( > 1 f or w > and f < 0, it follows that, 
if < p < 1 and £ < 0, this is also cm and E pQ (— z) > for 
(B3a) P > ^- Gi ven that -Cf^to[6(f)] = 1, we also find from this that 

(B5) 



£[d(t) + r 1 E-{(-tP)] = (l+w- p f. 

t—HV 1 ' 

For the specific discussion concerning sufficient conditions 
for the non-negativity of the df, we basically consider 

P(Y)R(r 2 ) = PQ¥)R (r z ) + r^PQV) 

where Rq(x) - R(x) - x~ l . 
would be 



r(8,L 2 ) = £- 1 



The corresponding df (£q = 0) 

Z&PiS). 



V27T1 



-S(Z/), 



and thus it is obvious that corresponding sufficient condi- 
tion is together + <9i P > and those derived in Sect. [6] with 
^o(^)- With R(x) given by equation (190b| i, the preceding dis- 
cussion in Sect. 17.21 actually extends to b = thanks to the 
non-negativity of E p Q (— z) > 0. It follows that all the corol- 
laries actually hold inclusively for /3\ = 1. [Note the con- 
dition +&t A P > and P(0) = ■ ■ ■ = p(W-D (0 ) = o for 

3 A > | -fa > | -ySi = | implies ^-P > 0.] 
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